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SUMMARY 


Nonlinear  beam  kinematics  are  developed  and  applied  to  the  dynamic  analysis  of 
a  pretwisted,  rotating  beam  element.  The  common  practice  of  assuming  moderate  rota¬ 
tions  caused  by  structural  deformation  in  geometric  nonlinear  analyses  of  rotating 
beams  has  been  abandoned  in  the  present  analysis.  The  kinematic  relations  that 
describe  the  orientation  of  the  cross  section  during  deformation  are  simplified  by 
systematically  ignoring  the  extensional  strain  compared  to  unity  in  those  relations. 
Open  cross  section  effects  such  as  warping  rigidity  and  dynamics  are  ignored,  but 
other  influences  of  warp  are  retained.  The  beam  cross  section  is  not  allowed  to 
deform  in  its  own  plane.  Various  means  of  implementation  are  discussed,  including  a 
finite  element  formulation.  Numerical  results  obtained  for  nonlinear  static  problems 
show  remarkable  agreement  with  experiment. 


1.  INTRODUCTION 


It  is  now  widely  recognized  that  aeroelastic.  analysis  of  helicopter  rotor  blades 
particularly  of  hingeless  and  bearingless  rotor  blades,  requires  the  incorporation  of 
kinematical  nonlinearity  (ref.  1).  The  main  reason  for  this  requirement  is  that  the 
stability  and  response  of  such  systems  depend  strongly,  in  some  cases,  on  the  cou¬ 
pling  between  bending  and  torsion  motion.  This  important  coupling  cannot  be  obtained 
accurately  without  consideration  of  the  kinematical  nonlinearities. 

A  brief  history  of  the  developments  of  nonlinear  equations  of  motion  for  rotat¬ 
ing  beams  prior  to  1974  is  given  in  reference  1.  These  developments  are  generally 
concerned  with  slender  beams,  with  the  effects  of  shear  deformation  ignored.  Since 
1974,  the  major  contributions  to  this  subject  have  been  those  of  Kaza  and  Kvaternik 
(ref.  2),  and  Rosen  and  Friedmann  (ref.  3).  The  equations  of  motion  developed  in 
references  1-3  are  very  similar.  In  references  1  and  3  an  ordering  scheme  is  used 
that  limits  the  kinematical  development  to  moderate  rotations.  In  reference  2  the 
nonlinearities  are  limited  to  the  second  degree  in  the  displacement  variables.  These 
two  different  methods  of  specifying  "moderate  rotations"  in  references  1-3  are  vir¬ 
tually  equivalent  for  this  problem. 

In  a  general-purpose  analysis,  a  single  set  of  equations  is  desirable  -  one  that 
is  valid  for  all  values  of  the  equation  parameters,  within  some  range.  When  moderate 
rotations  are  assumed,  situations  can  easily  arise  in  which  the  solution  violates  the 
assumption  of  moderate  rotations.  One  example  is  the  case  of  a  thin  beam  for  which 


the  ratio  of  bending  stiffnesses  is  small  compared  to  unity.  In  order  to  avoid  this 
problem,  certain  ad  hoc  modifications  to  the  equations  of  reference  3  were  intro¬ 
duced  and  were  necessary  in  order  to  produce  the  excellent  correlation  obtained  in 
reference  4  with  experimental  data  for  large  displacements  of  an  end-loaded  canti¬ 
lever.  For  example,  magnitudes  of  the  beam  bending  and  torsion  stiffnesses  had  to 
be  specified  before  the  equations  could  be  put  into  final  form  for  solution. 

Ideally,  the  magnitude  of  parameters  in  the  equations  for  general-purpose  analyses 
should  not  influence  the  equations  themselves.  Such  an  ideal  is  evidently  not  pres¬ 
ent  in  the  ordering  schemes  of  references  1  and  3  or  in  any  arbitrary  a  priori 
restriction  to  second-degree  nonlinearity,  as  in  reference  2. 

There  are  other  shortcomings  in  the  equations  in  references  1-3.  For  example, 
the  effects  of  pretwist  are  not  treated  rigorously.  An  improved  treatment  of  pre¬ 
twist  effects  is  presented  in  reference  5  for  a  simplified  problem  involving  only 
torsion  and  axial  displacement.  Additional  work  is  required  to  incorporate  those 
analysis  techniques  into  a  general,  nonlinear,  bending-torsion-extension  analysis. 

In  reference  6  exact  nonlinear  kinematical  relationships  are  developed  and  additional 
insight  is  presented  concerning  relationships  among  the  equations  of  references  1-3. 
Finally,  in  reference  7  it  is  shown  that  inconsistencies  are  virtually  unavoidable 
in  ordering  schemes  based  on  displacements  and  rotations  when  the  magnitude  of  the 
torsion  rigidity  is  small  compared  to  bending  stiffnesses. 

Because  of  the  problems  with  kinematical  limitations  in  the  above  approaches,  it 
seems  appropriate  to  model  the  kinematics  of  a  slender  beam  without  resorting  to  an 
ordering  scheme  on  rotations  or  to  arbitrary  restrictions  on  degree  of  nonlinearity 
allowed  in  expressions  involving  displacement.  This  method  would  avoid  some  of  the 
limitations  of  previous  analyses  and  circumvent  nonrigorous  modifications  such  as 
were  found  necessary  in  reference  4.  In  the  present  work,  the  development  in  refer¬ 
ence  6  serves  as  a  foundation  along  with  some  important  observations  from  references  8 
and  9.  Rather  than  develop  the  partial  differential  equations  of  motion  for  this 
problem,  the  objective  is  to  develop  a  statement  of  the  principle  of  virtual  work  for 
dynamic  analysis  of  a  rotating  beam  element  with  Euler-Bernoulli  kinematics.  This 
statement  will  serve  as  the  basis  of  a  Ritz-type  modal  model  for  an  entire  rotor  blade 
or  of  a  single,  finite  element  of  a  rotor  blade  without  having  to  consider  the  partial 
differential  equations  and  the  natural  boundary  conditions,  which  are  available  from 
the  statement,  if  needed. 

Although  the  results  presented  here  are  for  the  case  in  which  shear  deformation 
is  neglected,  the  necessary  relations  for  including  shear  deformation  are  included  as 
part  of  the  development.  Similarly,  with  the  appropriate  constitutive  law,  effects 
such  as  orthotropy  or  anisotropy  could  be  incorporated.  The  initial  curvature  of  the 
elastic  axis  and  effects  associated  with  open  cross  sections  could  also  be  incorpor¬ 
ated.  Tlie  detailed  development  of  these  topics  is  reserved  for  future  extensions  of 
the  present  analysis.  The  present  development  proceeds  as  follows:  The  beam  kine¬ 
matics  are  developed  in  section  2  in  such  a  way  that  the  necessity  for  an  ordering 
scheme  is  obviated.  This  section  includes  the  development  of  strain-displacement 
relations  and  direction  cosines  of  the  nominal  cross-sectional  plane.  In  section  3 
generalized  forces  caused  by  internal  loads  are  developed  from  a  consideration  of 
the  strain  energy.  Expressions  for  generalized  forces  caused  by  inertial  and  g’avi- 
tational  loads  are  developed  in  section  <4  based  on  the  work  done  by  these  loads 
through  a  virtual  displacement.  Generalized  forces  caused  by  a  general  set  of 
applied  distributed  loads  are  developed  in  section  5  similar  to  the  development  in 
section  4.  A  finite  element  implementation  is  discussed  in  section  6.  Finally, 


numerical  results  obtained  from  a  finite  element  calculation  of  nonlinear  static 
equilibrium  are  presented  in  section  7. 


2 .  KINEMATICAL  DEVELOPMENT 


In  this  section,  the  kinematics  for  the  beam  element  are  developed  starting  with 
the  rigid-body  motion  of  the  cross-section  plane  as  described  in  reference  6.  Cross- 
section  warp  is  then  superimposed  on  the  rigid-body  motion  to  obtain  the  final  dis¬ 
placement  field.  Next,  the  strain-displacement  relations  are  developed  from  the 
displacement  field.  The  extensional  strain  is  then  assumed  to  be  small  compared  to 
unity.  This  assumption  is  used  to  simplify  the  orientation  description  and  moment 
strains  considerably,  without  sacrificing  accuracy. 


2.1  Development  of  Displacement  Field 


First,  consider  a  straight-beam  element  with  the  associated  coordinate  systems 
shown  in  figure  1.  It  is  assumed  that  the  motion  of  the  frame  F  is  known  in  an 
inertial  frame  I.  A  set  of  dextral  axes  Xj^,  i  =  1,2,3  is  assumed  to  originate  at 
F*,  the  origin  of  F.  The  Xj-axis  lies  along  the  elastic  axis  of  the  beam  element. 
Each  unit  vector  b^,  i  =  1,2,3  is  parallel  to  the  corresponding  axis  xi.  The  beam 
is  assumed  to  be  pretwisted  so  that  the  local-beam  cross  section  is  rotated  by  an 
angle  9  about  the  Xj-axis  at  any  point  on  the  Xj-axis.  At  x^  =  0,  0  is  defined 
to  be  zero  so  that  the  major  axis  is  parallel  to  the  X2-axis  and  the  minor  axis  is 
parallel  to  the  Xj-axis.  Consider  an  arbitrary  material  point  in  the  beam  prior 
to  deformation,  denoted  by  Mj.  The  position  vector  of  Mq  with  respect  to  F*  is 

_  Xgbj  +  cos  9  -  sin  6)b^  +  sin  0  +  cos  0)b2  (1) 


Denote  the  same  material  point  in  the  deformed  beam  by  M.  The  position  vector  of 
M  with  respect  to  F*  is 


MF*  F  F  P  P 

R  =  X  b,  +  u.b.  +  5  b  +  YCxlACC  ,5  )b 
i-3  1-1  a-a  3' 


(2) 


where  repeated  Latin  indices  imply  summation  from  1  to  3  and  repeated  Greek  indices 
imply  summation  from  1  to  2 . 

The  axes  and  ^2  sre  along  principal  axes  for  the  local  cross  section  at  a 
point  P*  on  the  elastic  axis  and  remain  so  during  deformation.  The  unit  vectors 
b^  are  parallel  to  the  cross-section  principal  axes  at  P*,  which  is  located  at  the 
origin  of  a  dextral  system  P.  The  unit  vector  bg  is  defined  as  b^  x  b^  and  is 
thus  normal  to  the  cross  section.  The  displacements  uj^,  i  =  1,2,3,  are  along  b?, 
respectively.  A  warp  displacement  field  has  been  added  vectorially  to  the  rigid-body 
component  of  the  position  following  Wempner  (ref.  9).  The  warp  amplitude  is  V,  and 
X  is  the  cross-section  warp  function. 


The  basis  vectors  for  the  P  axes  are  denoted  in  equation  (1)  by 


where  Cij  is  a  3  x  3  matrix  whose  elements  may  be  specified  by  any  of  several  sets 
of  parameters  (ref.  10).  All  three-parameter  descriptions  of  the  rotation  have 
inherent  singularities.  Classical  Euler  orientation  angles  (called  body-two  orien¬ 
tation  angles  in  reference  10)  have  singularities  at  values  of  certain  angles  equal 
to  zero.  The  Tait-Bryan  orientation  angles  (called  body-three  orientation  angles  in 
reference  10)  have  all  singularities  at  90°.  This  fact  makes  them  more  amenable  to 
descriptions  of  rotation  caused  by  structural  deformation  than  Euler  angles,  since 
the  case  of  no  rotation  would  correspond  to  no  deformation  and  Euler  angles  are 
singular  at  that  condition  (ref.  6).  To  avoid  singularities  altogether,  at  least 
four  parameters  are  required  to  describe  the  rotation.  One  of  the  best  known 
descriptions  is  the  set  of  Euler  parameters.  Euler  parameters  enable  rotation  to 
be  described  with  four  scalar  quantities  upon  which  the  direction  cosine  matrix  ele¬ 
ments  depend  quadrat ically .  It  is  possible  to  eliminate  one  of  the  Euler  parameters 
algebraically  and  derive  the  Rodrigues  parameters.  There  is  a  singularity  at  180° 
of  rotation  in  any  direction,  but  this  is  rarely  a  problem  in  deformable  structures. 
The  direction  cosines  are  simple  ratios  of  quadratic  polynomials  in  the  Rodrigues 
parameters.  Furthermore,  the  inverse  operation  (i.e.,  given  the  direction  cosine 
matrix,  find  the  parameters)  is  trivial  compared  to  the  same  operation  with  orienta¬ 
tion  angles  (ref.  10). 

In  this  paper,  both  Tait-Bryan  orientation  angles  and  Rodrigues  parameters  are 
considered.  The  primary  development  is  executed  with  the  orientation  angles  because 
of  the  simplicity  of  the  result  and  the  familiarity  of  the  method.  Rotations  free 
of  singularity  up  to  90°  are  completely  acceptable  in  helicopter  blade  applications. 
For  applications  in  which  the  90°  restriction  is  unacceptable,  the  kinematical 
development  for  Rodrigues  parameters  is  given  in  the  appendix. 

The  direction  cosines  are  here  expressed  in  terms  of  Tait-Bryan  orientation 
angles  -j  i ,  i  =  1,2,3  defined  as  follows.  Let  basis  vectors  b?  be  introduced 
beginning  with  b^  aligned  with  b^.  Then  perform  sequential  rotations  Gi  about 
b?,  i  =  1,2,3  until  b?  is  aligned  with  the  principal  axes  of  the  deformed  beam 

i  =  1,2,3,  as  shown  in  figure  2.  Other  sequences  of  rotations  are  possible,  but 
here  9]^-02~63  is  used.  Direction  cosines  of  the  local,  def ormed-beam,  cross- 
section  principal  axes  P  with  respect  to  F  when  expressed  in  matrix  form  are 


where  cj^  =  cos  0^  and  Si  =  sin  The  displacement  field  is  now  completely  spe¬ 

cified  although,  as  pointed  out  in  reference  6,  two  of  the  three  angles  0^  and  02 
can  be  eliminated  if  shear  deformation  is  neglected. 


2.2  Development  of  Strain-Displacement  Relations 

Now  that  the  displacement  field  is  determined  in  terms  of  Uj[  and  G^j^,  it  is  a 
straightforward  matter  to  calculate  strain-displacement  relations.  There  are  several 
strain  measures  that  could  be  used  in  this  type  of  analysis.  Almansi  strain  was  used 
by  Hodges  and  Dowell  (ref.  1)  in  their  preliminary  development,  and  Green  strain  has 
been  used  by  almost  everyone  else.  For  infinitesimal  strains,  however,  differences 
among  the  common  definitions  of  strain  are  small  and  are  ignored  in  this 


4 


development.  The  intent  here  is  to  develop  a  set  of  strain-displacement  relations 
that  will  be  linear  in  elongations  and  shears,  but  unrestricted  in  rotations  up  to 
changes  in  orientation  where  singularities  are  encountered.  Thus,  it  is  immaterial 
whether  it  is  Green  strain  or  Almansi  strain  that  serves  as  the  starting  point  since 
the  relations  will  be  simplified  for  small  strains  anyway. 


M  MF" 

The  independent  variables  of  the  vectors  R  “  and  R  (namely,  >  €2’  ^3) 

constitute  a  nonorthogonal  curvilinear  coordinate  system  identical  to  the  one  used  in 
reference  5.  The  steps  followed  in  reference  5  to  derive  a  set  of  strain  displace¬ 
ment  relations  from  these  vectors  are  as  follows: 


1.  Obtain  the  covariant  base  vectors  for  the  undeformed  state  from  equation  (1) . 

2.  Obtain  the  covariant  metric  tensor  for  the  undeformed  state  from  step  1. 

3.  Obtain  the  contravar iant  metric  tensor  for  the  undeformed  state  from  step  2. 

4.  Obtain  the  relationship  between  a  local  Cartesian  coordinate  system  and 
step  1  from  step  3. 

5.  Obtain  the  covariant  base  vectors  for  the  deformed  state  from  equation  (2). 

6.  Obtain  the  covariant  metric  tensor  for  the  deformed  state  from  step  5. 

7.  Obtain  the  Green  strain  tensor  from  half  the  difference  between  step  6  and 
step  2. 

8.  Transform  the  Green  strain  tensor  from  step  7  to  the  local  Cartesian  coordi¬ 
nate  system  using  step  4. 

9.  Ignore  elongations  and  shears  with  respect  to  unity  in  all  strain  components 
(i.e.,  discard  all  squares  and  products  of  elongations  and  shears,  leaving  each 
strain  component  to  be  a  linear  combination  of  elongations  and  shears). 

The  details  of  the  algebra,  although  lengthy,  are  straightforward  and  are 
omitted  in  this  report.  The  engineering  strain  components  resulting  from  these 
operations  are 


31  ^2^  ^*^3 

■32  ~  ^32  ^  ^^2  ~  ®  ) 


^33  =  S3 

+  (52S  -  +  ^(^3  -  e')' 


The  strains  at  the  reference  axis,  referred  to  as  force  strains,  are  given  by 


^3i  ^ij  ^*^3  j  ^  3  i 


where  is  the  Kronecker  symbol.  The  curvature-like  quantities,  referred  to  as 

moment  strains,  are  given  by 
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1 


where,  in  matrix  form, 


K  .  =  (D  .  .  6 

1  iJ  J 


■^2^3  0 

“  =  -^2^3  Cj  0 


0  1 


Derivatives  of  the  warp  function  are  denoted  by  =  3X/3^qj,  A  restricted  warp 
amplitude  :  =  K3  -  0'  is  assumed  here  instead  of  the  more  general  formulation  in 
terms  of  :  explicitly  as  in  reference  5.  This  assumption  results  in  the  neglect 
of  transverse  shear  in  the  outer  fibers  of  the  beam.  It  is  interesting  to  note  that 
the  force  and  moment  strains  are  very  similar  to  those  of  Reissner  (ref.  8)  except 
that  the  ones  above  are  expressed  in  terms  of  Tait-Bryan  orientation  angles  instead 
of  Rodrigues  parameters  as  in  reference  8.  The  moment  strains  have  the  dimensions 
of  curvature,  but  differ  from  curvature  by  a  factor  of  s',  where  s  is  the  length 
coordinate  along  the  deformed-beam  elastic  axis  (ref.  6)  and  (  ) '  denotes  the 
derivative  with  respect  to  X3.  The  moment  strains  also  closely  resemble  those  of 
Wempner  (ref.  9)  developed  for  small  deformation  of  arbitrarily  curved  beams. 


It  is  possible  to  eliminate  and  02  from  the  analysis  if  shear  deformation 
is  ignored,  resulting  in  an  Euler-Bernoulli  beam  model.  For  this  case  the  vector 
tangent  to  the  beam  elastic  axis,  3R^^*/3slr  _n»  must  remain  normal  to  the  local,  beam 

cross  section  during  deformation.  Thus,  from  equation  (2) 

MF*  I  F  P 

3R  /3s|  ^  =  [dj^OXg/as)  +  3u^/3s]b^  =  bj  (9) 


By  virtue  of  equations  (3)  and  (9) 


3Xj  9u^ 
^ai  =  ^3i  ^ 


Since  C^j  is  orthonormal,  =  1.  If  all  derivatives  are  expressed  with  respect 

to  Xg  instead  of  s,  an  expression  for  s'  is  obtained: 


From  equations  (10)  and  (4) 


1'  =  [u]2  +  +  (1  +  up2jl/2 


=  u'/s' 


(12a) 


-s,c„  =  u'/s' 


(12b) 


-1^2  =  (1  +  up/s' 


(12c) 


The  angles  0^  and  0^  can  be  eliminated  from  with  the  following  relations 

obtained  from  equations  (12) 


6 


i.qii.it  ions  (12),  when  substituted  into  the  expressions  for  force  strains,  yields 


e,  =0 

3a 


‘'33  ® 


(14) 


Cons  idot  il>  1 V  more  algebra  is  required  to  eliminate  6^^  and  9^  from  the  moment  strains 
ii  thi  above  metliod  is  used.  it  is  relatively  simple,  however,  to  write  them  directly 
I  r  '111  t  eiefeiue  s  in  tet  ms  of  derivatives  with  respect  to  s  where  ()'*’=  3/9s(  ). 


(1 


++ 


-  n  .  ) 


(1  -  U+‘  -  U+2)1/2 

1  2 


(1 

(1  - 


b 

-  u+'O^/^s  u+u+u-^l 

_ 1_ _ 3 _  ++  1 . 

+2  +2-,  1/2  2  +2  I 

Ui  -  U2  )  I  1  -  J 


(15a) 


(15b) 


W1 

1  -  uf  _ 


(15c) 


i  .  o  .  '  i  I  !  'm  .  xi'rossed  in  terms  of  (  )'  quantities  instead  of  (  )"*"  quanti- 

t  .  .  ,11  M  i  n  (  i  I  )  iiul  leiiclhv  algebra.  liiis  is  unnecessary,  however,  in  light  of 

t:a  --i:',,  i  i  ;  i '  a  I  i ' '11 S  in  liie  next  section. 


J.i  Snial  1 -St  rain  S  impJ  if  icat  ion  of  Kinematics 

IVaanse  the  moment  strains  are'  complicated,  it  is  useful  to  simplify  them 
throiudi  the  derivation  of  a  sm.il  1-strain  approximation.  This  is  accomplished  by 
nec,Ieiting  the  1  onp,  i  t  ud  i  na  1  strain  ot  the  elastic  axis  with  respect  to  unity  in  the 
direction  cosines  as  well  as  the  moment  strains.  The  expression  for  the  longitudinal 
force  strain  s'  -  1  is  already  linear  in  the  elongation  of  the  elastic  axis;  the 
moment  strains  that  are  to  be  obtained  .are  then  independent  of  elongation  of  the 
elastic  axis  (i.e.,  independent  of  s').  It  is  interesting  to  note  that  Reissner's 
moment  strains  differ  from  curvatures  bv  a  factor  of  s'  and  are  independent  of  s'. 
When  u,  and  are  eliminated,  additional  s'  terms  are  introduced  into  the 


equations  through  equations  (13).  Simply  setting  s'  =  1  in  all  places  that  it 
occurs  in  the  moment  strains  and  direction  cosines  removec  the  dependence  of  these 
quantities  on  £33* 


Another  way  to  view  this  approximation  is  to  expand 
series  with  elongations  caused  by  stretching  of  the  elast 
eter.  This  separates  elongation  into  components  on  and  o 
only  the  terms  linear  in  elongations  are  retained,  the  re 
set  s'  =  1  in  all  quantities  except  £33-  Thus  strain 
to  unity  as  in  the  example  given  in  reference  6,  p.  32. 
that  this  sort  of  approximation  invalidates  the  strain  fo 
sional.  In  reality,  a  simplification  of  this  sort  cannot 
accuracy  of  the  mathematical  model  as  long  as  the  strains 
which  they  must  be  for  applications  of  Hooke's  law.  The 
the  derivation  and  in  the  final  equations  are  substantial 


the  strain  in  a  Taylor 
ic  axis  as  the  small  param- 
ff  the  elastic  axis.  When 
suit  is  equivalent  to  having 
is  considered  small  compared 
It  is  claimed  in  reference  11 
r  cases  other  than  inexten- 
significantly  affect  the 
are  small  relative  to  unity, 
resulting  simplifications  in 


is  set  equal  to  1  in  equations  (13)  and  (15),  the  result  is 


(1  -  ui^) 


(1  -  u'^)^/2c  r  u'u'u""l 

- - - -  u"  +  — 1—2—1^ 

1  -  ui^  -  ul^)^/2  2  1  _  ^,2j 

(1  -  u’2)i/2s  r  u'u'u"! 

- 1 - 3_^  L''  + 

1  -  u{^  -  u^^)^/^  ^  1  -  u^^J 


2x1/2 


(1  -  uj 

(1  - 


u '  r  u  'u  'u"  "I 

1  _  ,  "  +  12  1 

Ui^  -  uf^TI  ^*^2  ^ 


(16a) 


(16b) 


(16c) 


(1  - 


(17a) 


(1  - 


i2x 1/2 


(17b) 


(17c) 


=  (1  -  u’^) 


2x  1/2 


(17d) 


Note  that  the  sign  of  the  square  root  quantities  becomes  ambiguous  when  9^  or  0, 
exceeds  00°.  ITius,  it  is  imperative  in  this  formulation  to  restrict  orientation 
angles  te  less  than  90°.  Also  note  that  -  q  and  Cj j  are  now  independent  of  u'; 
s'  must  net  he  set  equal  to  unity  in  the  force  strain  Zix  =  s'  -  1. 


Geometric  boundary  conditions  are  determined  from  specified  u^  and  Cj^j  at 
Xj  =  0  and  =  •.  Natural  boundary  conditions  may  be  obtained  from  the  complete 

statement  ol  the  principle  of  virtual  work  as  developed  in  the  next  three  sections. 
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3.  DEVELOPMENT  OF  GENERALIZED  FORCES  CAUSED  BY  INTERNAL  LOADS  FROM  STRAIN  ENERGY 


In  this  section,  the  internal  loads  for  the  beam  element  are  developed  from  the 
strain-displacement  relations  from  section  2  and  a  conventional  constitutive  law. 

The  beam  is  assumed  to  be  of  such  a  configuration  that  open-cross-section  effects, 
such  as  warping  vljidit}j,  are  negligible,  which  is  justified  for  rotor-blade  cross 
sections.  However,  the  effects  of  warping  are  retained  in  other  parts  of  the  inter¬ 
nal  loads  development.  This  assumption  is  helpful  in  a  finite  element  context  since 
enforcing  kinematical  boundary  conditions  on  the  warp  displacement  field  at  finite 
element  nodes  is  not  possible  in  the  general  case  of  beams  being  joined  at  arbitrary 
ptusitions  and  orientations  with  respect  to  one  another. 

The  virtual  work  on  the  internal  loads  is  obtained  from  the  variation  of  the 
strain  energy  U  as  in  (ref.  9) 


5U  = 


( II 


where 


(Et,  fic,  +  Ge  6e  )dC  d?  dx 
3  3  3  3  3a  3a'^  ^2  3 


(18) 


3  3 


=  5s’  +  -  C  <5k  +  [(C  +  C  )(<  -  0’)  +  (5  X  -  C  X  )G']6k  +  X6k 

^112  123  2112  3 


>  (19) 


'5t  =  (X,  +  S.)5k 

i  c  ^  13 


J 


The  expressions  are  greatly  simplified  if  the  variations  of  force  and  moment  strains 
are  written  as 


as'  , 

os  =  T - T  Sul 

du^  1 


aK. 

5k  ■  =  •;;; — 'Su'.'  +  +  r — r  6u '  +  e.  K  66 

1  9u'  a  13  3  9u'  a  laa  a  3 

a  a 


where  is  the  Levi-Cevita  permutation  symbol  and 


9s' 

9u  . 
1 


aK 

a 


Ik 


dU  , 


■Si  +  "1 


a  2 


'3  3 


3k 

_ 3 

auy 


-c^  c 

3  1  32 


C  (1  -  ) 

3  3  3  1"^ 


3k  . 


-c . 

_ ^ 

c 

3  3 


u"  f  C-  ,C  C  C  I  u"C  C 

_ L  _  '  "  ‘  12  31  _  2  II  31  32 

cl  1  -  C'  ^33  (1  _  ' 

33*-  J  33'  ^31- 


(20) 


(21) 


(22) 


9 


'3k 

2 

u" 

1 

C 

32 

3u; 

L  - 

^  3  1 

1 1 

-u,  C 

11^;, 

(2c; 

c  c 

22  3  1 

Ss' 


u"C  C  C, 

2  21  31  32 

^33^^  -  Csx) 


Cj,)  u. (1  -  c: j 


3u; 


-  ^ 


C  (1  -  c'  )■ 


3k;  .  C  . 

1  11 


3u: 


(U1C31  +  u^Cj^) 


(22) 

(cont) 


'3  3 


Next,  the  strains  are  substituted,  equations  (5)  with  1 3  given  by  equation  (14), 
bv  equations  (16),  and  the  virtual  strains  given  in  equations  (19)  into  the 
strain  energy  bv  equations  (18).  The  resulting  expression  can  be  arranged  by 
terms  that  multiply  is',  and  SKj-  These  coefficients  are  denoted  by  F3, 

and  B,  respectively,  which  are  given  by 


F3  =  F(.(s'  -  1)  +  +  -y-  (K3  -  e’)^  +  DoO'(k3  -  9') 


t 


M.  =  i,.  .L.k  +  e.„,  F,(s'  -  1')  (k  -  0’)^  +  D. 

I  fS  C)  t  j  s  t  D  3  p  Z,  3  P 


5'<,:,  -  v')l 


B  3D3O' 

+  13(5'  -  1)  -  o_^^3B^K;g  +  2“  <^<3  - 


(ic  -  e’)  +  D  0’^ 

3  4 


] 


K  (K3  -  .')  +  [D^(s'  -  1)  -  e^^3D^K3lO- 


B  =  S^(.:3 


')  +  S 


(23) 


These  quantities  are  the  stress  resultants  given  in  terms  of  the  displacements.  The 
distributed  axial  fcjree  is  Fj;  local  distributed  moments  along  b?  are  denoted  by 
M|.  TIil:  distributed  binioment  is  B,  which  may  be  neglected  for  the  class  of  beams 
considered  herein.  Ihe  section  properties  used  in  equation  (23)  are  defined  by  the 
•‘ollowine  integrals  ovt..r  tlie  cross  section. 


=  JJ  E  dC^ 


>U  "  jj  I 

•'•'A  I 

)  E  if  EC  C  e  AC  dr  dC 

3  JJ^  ^Y’Y  '‘[-3  '1  ^2 


(24) 
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15.  Dowell,  D.  H. ;  Traybar,  J . ;  and  Hodges,  D.  H. :  An  Experimental-Theoretical 
Study  Nonlinear  Bending  and  Torsion  Deformations  of  a  Cantilever  Beam. 
J.  Sound  Vib.,  vol.  50,  no.  4,  Feb.  22,  1977,  pp.  533-544. 

lb.  Greenberg,  J.  M. :  .Airfoil  in  Sinusoidal  Motion  in  a  Pulsating  Stream.  NACA 
TN-1326,  Jan.  1947. 


TABLE  1.-  CANTILEVERED  BEAM  LOADED  WITH  TIP  WEIGHT 


Property 

Value 

Metric 

Standard 

E 

7.2919x10'°  N/m^ 

(10.576x10°  Ib/in.^) 

C 

3.022x10'°  N/m^ 

(4.383x10°  Ib/in.^) 

P 

2807  kg/m° 

(0.1014  Ib/in. °) 

li 

(Ec''t/12) 

l2 

(EctVl2) 

c 

1.270  cm 

(0.4999  in.) 

t 

.3178  cm 

(0.1251  in.) 

L 

50.76  cm 

(19.985  in.) 

TABLE  2.-  ROTOR  BLADE  PROPERTIES  (NACA  0012 
AIRFOIL  WITH  TWO  BLADES) 


Rotor  diameter,  m . 1.923 

Blade  length  (L)  ,  . 870 

Hub  offset,  %R . 9.51 

Chord,  cm . 8.64 

Taper . 0 

Twist . 0 

Maximum  tip  Reynolds  number  .  600,000 
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APPENDIX 


In  this  appendix  the  direction  cosines  and  moment  strains  are  expressed  in 
terms  of  Rodrigues  parameters.  To  remove  the  sign  ambiguities  of  the  development  in 
the  text,  the  direction  cosines  should  be  left  in  terms  of  u^,  i  =  1,2,3  so  that 


3  1 


=  (6 


3  1 


Up/s' 


(Al) 


This  is  not  necessary  in  the  text  because  the  use  of  orientation  angles  limits  the 
rotations  to  be  less  than  90°;  furthermore,  it  makes  the  computation  slightly  more 
involved.  For  the  use  of  Rodrigues  parameters,  we  make  use  of  similar  relationships 
derived  in  references  8  and  10.  The  direction  cosines  in  terms  of  Rodrigues 
parameter^  f .  are 

Cij  =  [(1  -  +  e..j^4>k]/(l  +  ^Jj^)  (A2) 

After  much  algebraic  manipulation,  can  be  eliminated  in  terms  of  the  third  row 

of  C 


<t>  C  )/(l 
3  30 


C  ) 

3  3 


(A3) 


The  use  of  equation  (Al)  then  yields 

“  (^e^gju^  +  4.3U(i)/(l  +  s’  +  up 


(A4) 


which  goes  to  infinity  only  when  the  beam  rotations  due  to  deformation  reach 
180°. 


The  moment  strains,  as  simplified  for  a  straight  beam  based  on  those  of  Reissner 
(ref .  8) ,  are 


which  can  be  expressed  in  terms  of  u^,  <j)  ,  u'p  and  (j)^  by  differentiation  and  sub¬ 
stitution  of  equation  (A4)  into  equation  (A5) .  It  should  be  noted  that  the  reference 
basis  used  by  Reissner  corresponds  to  the  principal  axes  of  the  local  undeformed 
beam,  instead  of  the  principal  axes  of  the  root  used  herein.  Equation  (A5)  thus 
differs  somewhat  from  Reissner's  equation  (49)  in  reference  8. 


remainder  of  the  structure.  Results  for  a  droop  angle  of  0°  and  precone  angles  of 
0°  and  5°  are  shown  in  figures  7  and  8.  The  agreement,  again,  is  very  good. 

CONCLUSIONS 


Nonlinear  beam  kinematics  for  small  strains  and  large  rotations  have  been 
developed  and  applied  to  the  dynamic  analysis  of  a  pretwisted,  rotating  beam  element. 
There  are  no  explicit  restrictions  on  rotation  caused  by  deformation  in  these  equa¬ 
tions—  only  the  extensional  strain  of  the  elastic  axis  is  required  to  be  small  rela¬ 
tive  to  unity.  The  only  restriction  on  the  magnitudes  of  the  orientation  angles 
used  in  describing  the  cross-section  orientation  is  that  they  remain  less  than  90°. 
For  applications  of  the  kinematics  where  larger  rotations  may  be  encountered,  a 
method  of  overcoming  the  restriction  on  the  magnitude  of  rotation,  which  utilizes 
Rodrigues  parameters,  is  presented  in  the  appendix. 

In  order  to  be  applicable  to  all  existing  rotor/hub  configurations  in  helicop¬ 
ters,  the  analysis  needs  to  be  incorporated  into  a  hybrid  multibody/f inite-element 
program.  This  incorporation  is  under  development.  Useful  future  extensions  include 
constitutive  equations  for  composite  beams  and  effects  of  shear  deformation,  warping 
rigidity,  and  initial  curvature. 


Ames  Research  Center 

National  Aeronautics  and  Space  Administration 

Moffett  Field,  California  94035,  January  18,  1985 
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To  obtain  the  discretized  matrix  substitute  the  shape  functions,  equation  (45),  into 
equation  (70)  and  integrate  over  the  element  using,  for  example,  Gauss-Legendre  quad¬ 
rature.  The  dimension  of  the  discretized  matrix  will  depend  on  the  number  of  internal 
degrees  of  freedom.  Its  contribution  to  the  system  in  terms  of  nodal  degrees  of  free¬ 
dom  is  straightforward  and  can  be  left  to  the  reader  to  determine.  Similarly,  ele¬ 
ments  of  the  gyroscopic  matrix  can  be  calculated. 


These  equations  can  be  programmed  either  for  finite-element  computation  as  they 
are  or  written  in  a  more  explicit  matrix  form.  When  the  equations  are  linearized 
about  static  equilibrium  determined  from  nonlinear  static  equations,  a  convenient 
approach  is  to  calculate  the  mass  and  gyroscopic  matrices  explicitly  as  above  and 
solve  for  the  stiffness  matrix  by  numerically  perturbing  the  total  static  general¬ 
ized  force. 


7 .  RESULTS 


In  this  section,  two  sets  of  numerical  results  are  presented  along  with  corre¬ 
sponding  experimental  data.  The  numerical  results  were  obtained  by  exercising  a  pre¬ 
liminary  version  of  a  multibody/f inite-element  program  presently  under  development  in 
which  the  beam  element  formulation  outlined  in  section  6  is  implemented.  The  examples 
were  set  up  for  calculation  using  only  the  properties  given  in  this  report. 

The  first  set  of  data  concerns  a  cantilevered  beam  loaded  with  a  tip  weight 
(ref.  15).  The  properties  are  tabulated  in  table  1,  and  the  experimental  configura¬ 
tion  is  described  in  detail  in  reference  15.  The  beam  was  modeled  with  one  element 
and  a  sulticiently  large  number  of  polynomials  to  achieve  convergence.  Results  for 
transverse  tip  deflections  u^  and  03  and  tip  rotation  63  are  shown  in  figures  4-6 
along  with  experimental  data.  The  agreement  is  good,  much  better  than  in  reference  15. 
This  good  agreement  is  achieved,  however,  without  ad  hoc  modifications  of  the  equa¬ 
tions  based  on  the  values  ot  beam  stil  1  nesses  as  is  done  in  reference  4. 

The  second  set  of  data  concerns  a  cantilevered  rotating  beam  in  axial  airflow. 

The  air  flowing  down  through  the  rotc^r  plane  is  induced  by  the  thrust  at  collective 
pitch  setting  y,  for  a  two-bladed  rotor.  Bending  and  torsion  moments  were  mea¬ 
sured  near  tlie  beam  root  for  variable  tlirust  settings.  The  program  was  set  up  to 
caKulate  air  heads  based  (ui  .1  quas  i-ste.idv  aerodvnamic  formulation  by  Greenberg 
(ref.  16)  and  on  uniform  induced  inflow  di*termined  from  momentum  theory.  (The 
experimental  configuration  is  described  in  a  forthcoming  publication  by  Sharpe 
(Sharpe,  I).  I..:  An  Experimental  Investigation  of  tlie  Flap-Lag-Torsion  Aeroelastic 
Stibilit'.-  of  .1  Sma  1  1 -Sea  1 1’  Hingeless  Helicopter  Rotor  in  Hover.  NASA  TP,  tt)  be  pub¬ 
lished  in  1185).  The  rotc'r  properties  are  t.ibulated  in  tables  1  and  1.)  The  rotor 
blade  was  modeUal  with  two  uniform  ihastic  si-gments  and  rigid  masses  ft'r  the 
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riiL'  virtual  work  of  forces  and  moments  at  the  root  and  tip  nodes  must  equal  the  vir- 
tu<il  work  caused  by  generalized  forces  at  the  root  and  tip  of  the  beam 
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where  R^  is  the  transpose  of  R  and  the  Q's  are  the  generalized  forces  from  the 
discretization  of  the  beam  (i.e.,  coefficients  of  the  6q's). 

The  element  mass  matrix  may  be  written  explicitly  from  equation  (39)  by  collect¬ 
ing  coefficients  of  ^u  and  After  simplification,  and  ignoring  rows  and  columns 
associated  with  frame  degrees  of  freedom,  the  mass  matrix  can  be  obtained  from 


CT'R 
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The  forces  and  moments  at  the  root  and  tip  nodes  must  be  tie  t  e  rm  i  ned  i  roni  tin 
generalized  forces  at  the  beam  root  and  tip,  respectively.  First  note  tb.it  t  lu  v 
tual  displacement  and  virtual  rotation  of  the  root  node  must  be  ideiuital  t.>  tbii^, 
of  the  beam  root 


6u  (t)  =  'tu(0,t) 


^''(t)  =  ^(0,t) 


where  ^  =  61;  (see  eq .  36).  Similarly  at  the  tip 


cu^(t)  =  6u(«,,t) 


£'/(t)  =  6^.(Ji,t) 


The  relations  for  the  virtual  displacements  are  simple  first  variations  of  equa¬ 
tions  (55)  and  (58).  The  rotations  are  more  involved.  The  use  of  equation  (36) 
evaluated  at  the  root  yields 


L  a 


+  6  .  5q 

3  1  ^  “I : 


C.  .(0,t) 


where  =  ocS-b?.  In  matrix  form 


=  R(0,t)  J  6ij^ 


where 


RlXj  ,  t ) 


1  -  cli 


A  similar  relation  holds  for  the  beam  tip 


Similarly,  the  rotations  must  be  the  same The  direction  cosines  of  R  with 
respect  to  its  undeformed  position  R,  are 

C^'^(t)  =  C(0,t) 

At  the  tip,  the  displacements  of  the  node  and  beam  are  identical 

T  T  F  F  F 

u  b.  =  q  b  +  q  b  +  q  b 
Tl-l  ^13-1  ^23-2  ^32-3 


and  the  rotations  of  the  tip  node  ^  are  identical  to  those  of  the  beam  tip 
C(il,t) 

T  *  R  T  TR 

C  =  C  (t)C  =  C(2.,t)  (59) 

Equations  (56)  and  (59)  each  represent  nine  equations,  but  only  three  are  independent 
Three  preferred  elements  to  equate  are  C2i>  and  C32  since 


C2j^  =  -(1  -  Uj^)^^^  sin  0 


C31  = 


^32  *^2 


and  are  thus  easily  expressed  in  terms  of  q's  at  the  root  and  tip.  The  use  of 
equations  (46),  (47),  and  (60)  yields  for  the  root,  from  equation  (56) 


*^02  ^  ^R'R 


q  =  -sin"^ 
^41 


(i  -  c?;*') 


2\l/2 


and  for  the  tip,  from  equation  (59) 


and  on  the  first  derivatives 


u'(0,t)  =  q  „(t)/Jl  ;  u'((!-,t)  =  q  ,(t)/£ 

a  '  ^a2 '  > '  •>  a  ^a4 

The  standard  linear  displacement  functions  satisfy  these  conditions  for 

=  1  -  X 

and  the  standard  cubic  displacement  functions  satisfy  these  conditions  for  T : 


(47) 


1  -  3x^  +  2x 


(49) 


^2  =  X  -  2x^  +  2x'‘ 

Tj  =  3x^  -  2x^ 

T.  =  -x^  +  x^ 

4 

If  exceeds  3  or  Nj.N^  exceed  1,  the  following  higher-order  shape  functions 

allow  extra  generalized  coordinates  to  be  introduced  (ref.  13)  but  still  to  fulfill 
the  above  end  conditions 


■i 


n+3 


=  x(l  -  x)Gr,(5,3,x) 


n  =  0,1, 


(50) 


Vn+5  =  -  x)^G„(9,5,x) 

where  G^  is  a  Jacobi  polynomial  (ref.  14).  At  the  ends,  the  displacement  and  rota- 


1  to 

generalized 

coordinates  from  simple  kinematics. 

Let 

R 

R  .  R 

T  T  ,  T 

u  = 

'  u„  .  b  .  ; 
Ri-i  ’ 

u  =  u„  .b  . 

Ti-i 

(51) 

,  R 

,  F 

,  T  TR,  R 

b;  = 
-1 

^  ^i  ’ 

b.  =  C.  b 
-1  iJ-J 

(52) 

TR 

where  C  depends  on  the  pretwist  0(11)  of  the  beam  element  at  the  tip  of  the 
e 1 ement . 


,TR 


"0  ° 


-"0  ^0  ^ 


0 


0  1 


(53) 


where  s..  =  sin  y(?)  and  Cy  =  cos  'HZ).  At  the  root  of  the  element,  the  displace¬ 
ments  of  the  node  and  beam  are  identical 


bt‘  discretized  in  several  ways.  One  way  is  to  use  a  set  of  admissible  functions  in 
a  Ritz-type  analysis  based  on  equation  (44)  where  all  the  frame  variables  are  pre¬ 
scribed  and  rI'*!  and  are  zero.  A  variation  on  that  type  of  development  would 

be  to  allow  *  1  and  ' ,  *■  ^  to  be  nonzero,  thus  collecting  frame  (i.e.,  rigid  body) 

forces  and  moments  from  equations  (39)  and  (43)  directly.  In  this  case  there  would 
be  prescribed  and  nonprescr ibed  components  of  each  of  the  kinematical  quantities  in 
equation  (36).  There  are  clearly  other  ways  of  using  the  frame  motion  to  advantage. 

If  a  finite-element  discretization  is  used,  it  should  be  based  on  variable-order 
shape  functions  (refs.  12  and  13)  to  avoid  poorly  approximating  the  geometric  stiff¬ 
ening  term  which,  in  this  anal  sis,  is  calculated  from  the  longitudinal  strain 
s'  -  1.  If  Uj  is  crudely  approximated,  this  term  will  be  inaccurate.  It  should 
be  noted  that  in  a  redundant  structure  undergoing  finite  deformation,  the  geometric 
stiffening  effect  rrust  be  calculated  from  the  strain-displacement  relations. 


In  a  finite-element  implementation  the  frame  motion  and  forces  may  be  used  to 
advantage  when  coupling  elements  together  that  are  defined  in  different  moving  coor¬ 
dinate  systems,  such  as  at  the  interface  between  rotating  and  nonrotating  components 
of  a  helicopter  or  at  a  hinge. 

In  figure  3,  the  beam  element  is  shown  with  the  frame  F  and  two  nodes  R  and  T 
at  the  root  and  tip  of  the  beam,  respectively.  The  displacement  and  orientation  of 


the  beam  cross  section  as  a  function  of 


is  represented  by  the  displacement  and 


rotation  of  the  nodes  R  and  T  and  by  a  variable  number  of  generalized  coordinates 
that  are  kinematically  uncoupled  from  nodal  translations  and  rotations.  The  beam 
displacements  at  the  root  and  tip,  then,  must  be  determined  from  the  nodal 
d isplacements . 

The  bending  displacements  u^  and  u^  ate  to  be  expanded  in  C^-type  functions 
Tj^,  and  Uj  and  63  in  C°-type  functions  6i,  namely 


u  =  V  q  .(t)T.(x) 
a  ai  1 


^3  =  L  q3i(t)6.(x) 


S  =  E  qm(t)B.(x) 


where  x  =  0  <  x  <  1,  and  and  are  the  orders  of  the  shape  function 

polynomials,  N.^  >  3  and  N'a.Ni,  >  1.  The  generalized  coordinates  qj^j^  with 


j  -  1,2,3  have  the  dimensions  of  length  and  q^ 
conditions  on  the  functions  are 


is  dimensionless.  The  boundary 


u^(0,t)  =  q_.(t)  ;  u^(£,t)  =  q(^3(t) 
u^(o,t)  =  93^(6)  ;  u^(s,t)  =  qj^Ct) 


:  (0,t)  =  q.  _ (t)  ; 


(?,,t)  =  q  (t) 


where  the  section  integrals  are  defined  as 


and 


m  =  jj  p  d^^  dC^  ;  ^3  ~  jj  ^ 

iHf  H  jy  pC  dC^  dC2  ;  i  ^  JJ*  *^^^2 


F.  .  ,  F  F..  ..  ,F 

u=u.b.  ;  u=u.b, 

1-1  -  1-1 


5.  DISTRIBUTED  APPLIED  LOADS 


(4C 


(4] 


In  this  section,  generalized  forces  caused  by  distributed  applied  loads  are 
simply  stated  for  completeness.  The  applied  force  F  and  mo’^ent  M  act  on  the 
elastic  axis  so  that  the  virtual  work  is 


-6W 


app 1 ied 


loads 


r  ^  F*  I  F 

-I  {F  •  +  '^u)  +  M 

*0 

+  [(Xj  +  u)  X  F]  •  ^^^}dXj 


(4; 


or 


-6W 


applied 


loads 


Ok. 

'Sui  + 

p  F 

lb  dx,  -  F  •  (Su 

/' 

1  ^^3 

•  j  [M  +  (x-  +  u)  X 

F]dX3 

(4 

•'o 

•'o  "  " 

As  witli  the  inertial  loads,  the  generalized  forces  can  be  obtained  from  the  virtual 
work  and  anv  suitable  discretization  of  u-  and  0,. 

1  v 


6.  VARIOUS  SCHEMES  OF  IMPLEMENTATION 


The  total  virtual  work  may  be  written  as 

6U  -  6W.  ,  -  6W  .  .  ,  ^  =  0  (4- 

body  applied  loads 

where  the  various  components  of  the  total  virtual  work  may  be  obtained  from  equa¬ 
tions  (25),  (30),  and  (43).  The  kinematical  variables  u,-,  i  =  1,2,3,  and  0,  may 


■  (S; 


The  angular  acceleration  of  P  in  F,  ,  is  simply  the  time  derivative  in  F  o 
,PF  _ 


PF  F.PF 
a  =0) 


/ 3k  .  3^k  .  \  p 

=  (^  -a  +  ^ai^'a  + 


The  generalized  forces  are  now  found  from  the  virtual  work 
-“body  =  /'  //  “<5""  -  l>  •  “3 


'0  "A 

J  JJ  •  '  ~  +  t(23  +  H  +  §) 


MT  FT  MT  PF 

X  (a  -  g)]  •  ^  +  [5  X  (a  -  g) ]  •  _^  }dc^  d^^  dx, 


The  generalized  forces  associated  with  frame  motion  are  useful  in  multibody/finite 
element  applications  in  which  it  is  necessary  to  couple  bodies  together  that  are 
moving  relative  to  each  other.  These  generalized  forces  are  the  coefficients  of 
'Rp£  and  O'lpp  and  can  be  obtained  in  a  straightforward  manner;  the  calculation  of 
these  forces  is  left  to  the  reader.  For  the  generalized  forces  associated  with 
'Uj^,  and  the  necessary  ingredients  are  suitable  discretizations  for  Up 

and  Cj  and  the  virtual  work,  which  is  given  by 


■“body  -  !: 


F  r  ^ 

b.  •  I  6u. 

-1  J.  1 


+  a  ^  X  (x^  +  u)  +  X  X  (x^  ■+■  u)  ]  + 


F-  F..,  FI  FI  FI  FI  PF  FI  PF 

X  u  +  u)  +  a  X  m.  +  ^  X  (w  x  m^)  +  2w  *  5},- ~  2to  •  w  m. 

^  PF  ^  ^  PF  .  PF  ^  PF  FI  PF  FI  \ 

+  a  ^  E},-  -f  j  X  (0;  X  m,)  +  m  •  m_'jj  -  w  •  co  m^Ndx^ 


i  ■  ^  +  '^3i'^®3)(Hl,-  ^  -  g  +  X  (X3  +  u)  -f 

'3^^  '<  (X,  -f  u)]  +  X  +  i.u^^  -  •  1  +  •  i)  X 


We  must  find  the  acceleration  and  virtual  displacement  of  the  point  M  where  the 
position  of  M  with  respect  to  F*  is  governed  by 


MF* 

R  =  X3  +  u  +  5 


(27) 


where 


.  F 

^3  "  ^3^3 


u  =  u  .b  .  > 

1-1  ' 


(28) 


5  = 


a- a 


J 


The  velocity  of  M  in  an  inertial  frame  is  then 


MI  F*I  ,  FI  ,  ^  ,  F.  PI  ^  ,  P; 

V  =  V  4-  w  X  (x^  +u)+  u  +  o)  5 


(29) 


T7  •  * 

where  “^  (  )  is  a  time  derivative  in  F.  Cross-section  distortion  5  is  to  be 
neglected.  The  angular  velocity  can  be  written  as 


PI  PF  ,  FI 
w  =  u  +  w 


(30) 


Differentiation  yields  the  angular  acceleration 


PI  FI  .FI  PF  ^  PF 

a  =a  +w  xo)  +ct 


(31) 


Now  the  acceleration  can  be  written  in  terms  of  previously  developed  expressions  as 


MI  F*I  FI  FI  FT  .  FT  F  F  PT 

a  =  a  +  a  x  (x^  +  u)  +  w  x  [u  x  (x^  +u)]+2w  x  u+  u+a 


u.  PI 

+  a  X  I  4 


"-3 

PI 


/  PI 

X  (  ..,  X  t.) 


(32) 


Substitution  of  equations  (30)  and  (31)  into  (32)  yields 


MI  F*I  .  FI  ,  ,  ,  ,,  ,  FI  .FI  /  ,  ,  rM  ,  o  FI  F.  ,  F- 

a  =  a  4  a  x  (x^  4u45y+w  x[u  x  (x^  4u4£)]4  2w  x  u4  u 


FI  PF  PF  PF  PF  FI  PF 

4  2w  X  X  5)  +  u)  X  (u  X  O  +  w  X  (o)  X  5)  4  a  X  5 


(33) 


The  virtual  displacement  is  obtained  by  the  replacement  of  ^(’)  with  1^6 (  ),  a  vari¬ 
ation  in  the  frame  F,  in  equation  (29)  and  the  substitution  of  6R  and  6i|/  for 
V  and  respectively,  in  equations  (29)  and  (30),  as  in  reference  6,  which  yields 


Ml  F*I  FI  F  PF 


(34) 


I'*  I  FIF  FT  FIF 

where  ^R  =  (5Rplb£  and  ^  =  6i^p(bf  .  These  quantities  reflect  any  nonprescr ibed 

motion  of  the  frame  in  inertial  space.  The  angular  velocity  and  virtual  rotation  of 


the  cross-section  axes  P  in  F  (ref.  6)  can  be  written  as 


>1  CD 


I2  =  JJ^  dC,  d?. 


D4  =  JJ  E(e,33X^5g)^  dC^  d^ 


-  jj  E5»«b5b  ■‘Si  ‘‘S;  1  ■>  -  ‘^'<'1  +  *2>“  *  <^2  '  >,)"l<i5i  K, 

»3  -  ff  ■=<'•««»)'  <*«3  '*«2  •  Jf  ■‘f-i  'le,  =  »  ;  Ij  = 

A  A 

S2  =  EA(C,X^  -  ?,x^)d^  d5, 

S.  =  JJ^  EXCa^cc  dCi  d?2 

JJ  EXC^  d5^  d^  =  0  ;  JJ  EX  de^  d^  =  0 


(24) 

(cont) 


The  strain  energy  is  then 


^  6u^  +  M36e’jdx3 


+  M.  6u^  +  M36eydx3  (2f 

A  suitable  discretization  of  and  63  is  sufficient  to  obtain  generalized  forces 
in  a  Ritz-type  formulation. 


4 .  GENERALIZED  FORCES  CAUSED  BY  INERTIAL  AND  GRAVITATIONAL  LOADS 


In  this  section  the  generalized  forces  caused  by  inertial  and  gravitational 
loads  (i.e.,  body  forces)  are  developed  for  the  beam  element.  The  effects  of  cross- 
section  warp  on  the  body  forces  are  negligible  and  are  not  considered.  It  is 
assumed  that  the  gravity  field  g,  the  velocity  of  point  F*  in  an  inertial  frame 
E*^,  the  acceleration  of  point  F*  in  an  inertial  frame  a^*! ,  the  angvilar  velocity 
nd  angular  acceleration  of  frame  F  in  an  inertial  frame  and  ,  respectively 

are  given  and  defined  in  frame  F  as 
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